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Abstract. 

In this paper we investigate the question how spontaneous symmetry breaking 
works in the framework of Loop Quantum Gravity and we compare it to the results 
obtained in the case of the Proca field, where we were able to quantize the theory 
in Loop Quantum Gravity without introducing a Higgs field. We obtained that the 
Hamiltonian of the two systems are very similar, the only difference is an extra scalar 
field in the case of spontaneous symmetry breaking. This field can be identified as the 
field that carries the mass of the vector field. In the quantum regime this becomes a 
well defined operator, which turns out to be a self adjoint operator with continuous 
spectrum. To calculate the spectrum we used a new representation in the case of the 
scalar fields, which in addition enabled us to rewrite the constraint equations to a finite 
system of linear partial differential equations. This made it possible to solve part of 
the constraints explicitly. 
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1. Introduction 

Currently spontaneous symmetry breaking is the most accepted tool to define mass 
to particles. Its success can be observed especially in the case of vector fields since 
their original Lagrangian - the Proca Lagrangian - is non-renormalizable. In [lj we 
showed how one can quantize the massive vector field in Loop Quantum Gravity 
without spontaneous symmetry breaking. The main problem was that the Proca field 
had a second class constraint algebra which made it almost impossible to apply the 
framework of LQG. But with the help of symplectycal embedding one could eliminate 
these difficulties. Now the question arises what is the difference between the two 
theories. To study this one first has to apply the framework of LQG to a system where 
spontaneous symmetry breaking is used to generate mass for a U(l) vector field. This 
is done in sections 2 - classical theory and 3+1 decomposition - and 3 - regularisation 
and quantization. In section 4 we introduce a new basis for the scalar fields which is 
motivated by the fact that these are eigenstates of the configuration variables. It turns 
out that with the help of this new we are able to (partially) solve the constraints of 
the theory - this is done in section 5. Here we also analyze special solutions in order 
to understand the role of the scalar field. In particular we find that some of these 
are almost identical to the solutions obtained for the Proca field, thus we are able to 
relate the two theories. Section 6 deals with the "mass operator" and its properties, 
concentrating especially on those cases where the eigenvalues of this operator can be 
identified with the mass parameter of the Proca field. We obtain that in a sense the 
case of the symmetry breaking is a linear combination of infinite Proca theories. 

2. Classical theory 

In this section we will analyze the general framework of spontaneous symmetry breaking 
from a Hamiltonian perspective. In the first subsection we derive the 3+1 decomposition 
of the theory, while in the second we first have the U(l) field a VEV then perform the 
3+1 decomposition (this is useful because the similarities between the Proca field and 
spontaneous symmetry breaking become more transparent). 

2.1. Symmetric theory 

For simplicity we use the Lagrangian of a U(l) vector field (electromagnetic field) 
coupled to gravity and a U(l) complex scalar field on a space-time manifold M. The 
Lagrangian of the matter part is 

jjnat = I # x J=jj 
J M 

where 

Z)( 4 )$ = (8^ +ieA^)§ 



(-^ife (4)afe " \D^^D^ - \^ - a 2 ) 2 j , (1) 
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and (j, and a are positive constants. To distinguish between the electromagnetic field 

and the gravitational field the variables of the former are underlined. 

First let us make the 3+1 decomposition. Introduce on the space-time manifold M a 

smooth function t whose gradient is nowhere vanishing and a vector field t a with affine 

parameter t satisfying t a V a £ = 1. This gives a foliation of space-time, i. e. each t defines 

a 3-dimensional hypersurfice E t . Let us decompose t a into its normal and tangential 

part 

t a = Nn a + N a , (2) 

where n a is the unit normal of the hypersurfice H t , iV is the lapse function, N a is the 
shift vector. Define the induced, positive-definite metric on H t via 

q a b = g a b + n a n b (3) 

As it was done in [2] and pQ we define the pull-backs = q^A^ , D a = q b D^ and define 
Ao = t a A^, = t a A^. Substituting these into the Lagrangian one obtains 

L™< = J dtj d?x 2" ) - jJV9f(*** - « 2 ) 2 - 



,,, n n rU - ieA°$* - N a (D a $)*)(£ t § + ieA°$ - N b D b $) 



q ca D c $*D d <5> - 



N 2 

(4) 

where 

Ea = } § (£4, - D cA° ~ e abc B b N c ) (5) 
is the electric field and B_ a is the magnetic field. We now define the canonical momenta 

5L A$* - ieA°<!>* - N a (D a <S>Y 

* =JZ^ = ^ — (7) 

6L C t <$> + ie£<$>-N b D b <$> 

71 = Jc^ = ^ q N • (8) 

Finally we perform the Legendre transformation to arrive to the Hamiltonian: 

H mat = J d 3 x(NH mat + N a H™ at + A G) (9) 

n mat = q ab ( -- 2 ^_-- ) + + \y/q{V a *yV a * + ~ a 2 ) 2 (10) 

HT l = e abc B c E b + 7rX? a $ + 7r*(X> a $)* (11) 

G = DaEa + 2e(7r*$* - 7T$) 

where 7i ma * , 7i a aat and G are the matter contributions to the Hamiltonian- and 
diffeomorphism (modulo gauge transformations) constraints, and the electromagnetic 
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Gauss constraint. 

The (non-smeared, non-trivial) Poisson-brackets are: 

{tt,$} =8{x,y) (12) 
=6(x,y) (13) 
{E a ,A b } = 6 b J(x,y) (14) 

Since we are going to do symmetry breaking with the help of the scalar field, we write 
here the transformation rules for the scalar fields and their canonical momenta with 
respect to infinitesimal gauge transformation: 

{G(A),$} = -ieA$ 
{G(A),$*} = ieA$* 
{G(A),7r} = ieA% 
{G(A),tt*} = -ieAn* 

Before we continue, there are a few interesting observations that should be 
mentioned here: 

- All the constraints are real and only the scalar fields and their canonical momenta 
are represented by complex variables (note that in the Hamiltonian picture $ and 
$* are independent variables). 

- The transformation $ <-> $*, n <-> 7r* is a canonical transformation. 

- The true diffeomorphism constraint 7i™ ai + A^G is independent of the coupling 
constant e (it contains partial derivatives only). 

- This system has a first class constraint algebra, further more all the components of 
Hmat are gauge invariant respectively. 

2.2. New variables 

In spontaneous symmetry breaking first we introduce new fields r\ and in the following 

way: 

$(x) := (a + t)(x)) exp (i^) • (15) 

These variables are useful because the U(l) symmetry of the theory becomes more 
transparent. If we substitute this into the Lagrangian, we obtain 



J mat 



= J d*x^-g \-\f<$fW - i^V 4 N- 

(16) 



( a + v f __ + eA a )[ — + eA a \- - m \2a + V ) 2 



\ a \ a / 4 



If we compare this with the action of the symplectically embedded Proca field we 
immediately recognize the similarities between the two theories. The main difference is 
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that where the Proca theory had a parameter (m), now we have a field (rj + a). We 
wish to see how the Hamiltonian looks like in terms of the new variables. To do this, 
first we do the 3+1 decomposition of the above Lagrangian. Repeating the steps of the 
previous section first we define the canonical momenta 

s ° = ^£ = ^ (17) 

- SL (18) 



SCtV N 



5L f a + r,\ ( £ t Q - N a d a G + aeA - aeN a A a 



sc t e V a J ^ ) (19) 

and after the Legendre-transformation we obtain the Hamiltonian 

H mat = J d 3 x(NH mat + N a H™ at + A G) (20) 
E a E b + B a B b \ irl 1 



n mat = g ab ( — + ^daVdaV 4 



+ 



a + V J 2^ + ^(^) (9 a Q + aeA a f + y q ^2a + nf (21) 

H™ at = e ahc B c E h + -K v d a r] + n e d a O + aeA a ii e (22) 
G = - ae-Ke (23) 

The (non-smeared, non-trivial) Poisson-brackets: 

{*v,v} =\Kw) ( 24 ) 
{vr e ,e} = \8{x,y) (25) 
{E a ,A b } = 5 b J(x,y) (26) 

Now if one compares the Hamiltonian (1201) with the original (Q, it is easy to see 
that the two are connected with the help of the following canonical transformation: 

<P := [a + rj) exp I — 



a 

iG 
a 

ia \ f iQ 

- — vr© exp 

a + r] J \ a 

ia \ f iQ 



$* := (a + rj) exp 

/ ia 
7T := I 7r v - — — 7i e J exp ( - 

it* := ( ix T H ■ — 7r e I exp I — 

V a + rj J \ a 



Spontaneous symmetry breaking in Loop Quantum Gravity 



6 



Further more the above system is very similar to the case of the symplectically 
embedded Proca-field. To see this, let us introduce the canonical transformation 
7T0 — > — , — > eaO and define m 2 = e 2 (a + rf) 2 . Then we will obtain exactly 
the Hamiltonian of pQ, with the exception of a potential term. There are two major 
differences: there is an extra dynamical scalar field in the theory and the "mass" is 
constructed from the field r\. The later will be quite important since after quantization 
all the fields will become operators so we can define a "mass operator" , which spectrum 
can be identified as the mass spectrum (in pQ the mass was a parameter of the theory). 

2.3. Classical symmetry breaking 

In quantum field theory we use the unitary gauge to do gauge fixing. In the U(l) case 
this means we introduce the gauge-fixed vector field 

£\x):=A^\x)-±d^Q(x). (27) 
Substituting this into the Lagrangian we get 



Lmat 
1 



d A X\/—g 
fir] 2 (2a + ri) 2 



(28) 



Again we want to see how the Hamiltonian changes, so we do the 3+1 decomposition 
as we did in the previous sections. The canonical momenta will be 



7tr. 



5L 



SC t A 



E a 



v ~ 5C tV 
and the constraints will be 



5L Ct V - N a d aV 

V9- 



E a E b + B a B b , n 2 1 

H mat = q ab — + — + -VqdaVdaV + 

i i a _ i\ja a 

+ \^ V 2 (2a + rj) 2 + \e 2 ^{a + V ) 2 (A a A a + C ^ f) (29) 



HT at = e abc B c E b + Tc r d a Tj + e 2 A a ^(a + rj) 
G = 



2 Ap-N-Aq 

N 



9 ,-, ^Ao-N a A„ 
D a E a -e 2 ^(a + r ] ) 2 - 



N 



(30) 
(31) 
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R,t?} =^S(x,y) (32) 
{E a ,A b } = 5 b a 5(x,y) (33) 

If we compare this gauge fixed Hamiltonian with (120]) . we can see that in the 
Hamiltonian formalism the gauge fixing is equivalent to the introduction of the following 
two constraints 



C a := d a Arg($) = (34) 

C := a/Kf - e^-q(a + v) 2 ^'^ = 0, (35) 

(the second is equivalent to £ 4 6 = 0). This is precisely the gauge we used in the 
case of the symplectically embedded Proca field. There we showed that in LQG gauge 
fixing is not necessary, in fact it makes the quantization extremely difficult if not impos- 
sible. So we will not fix the gauge, instead we will try to solve the constraint related to it. 

To conclude we summarize the most important results of this section. 

We checked how one can implement spontaneous symmetry breaking in the Hamiltonian 
formalism. It turned out that introducing new variables means a canonical transforma- 
tion, while gauge fixing (as it was shown earlier e.g. in [13]) can be done by introducing 
new constraints. Interestingly these are exactly the same conditions which were intro- 
duced in the case of the Proca field in p]. Further more the Hamiltonian (|20|) is very 
similar to the symplectically embedded Proca Hamiltonian ([1], page 5), the only two 
exception is that we have an extra scalar field and the mass is not a parameter but 
defined with the help of the field rj. 



3. Quantization 

3.1. Gauge fields 

Quantization of the gravitational and electromagnetic field can be treated on the same 
footing since both are gauge fields - the gauge group of the former is, in the Ashtekar 
variables (0), SU(2) while the latter is a U(l) field. The detailed analysis of the method 
can be found in [3]-[8], here we just sketch the main idea and the notations. 
Let us consider a Yang-Mills gauge field with a compact gauge group G. The Hilbert- 
space can be constructed in the following way: let 7 be an oriented graph in X with 
ex,...,SE edges and V\,...,Vy vertices. Let h ei be the holonomy of the G- valued 
connection of the field evaluated along the edge. Let us define a cylindrical function 
with respect to a 7 graph in the following way: 
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f y (A):=f(h ei ,...,h eE ) (36) 

where / 7 is a complex valued function mapping from G E . The Hilbert-space of 
the Yang-Mills field is defined as the set of all cylindrical functions which are square- 
integrable with respect to a suitable measure (the Ashtekar-Lewandowski measure): 

H:= L 2 {A,dfi ALtG ) (37) 

In our case, G = SU(2) x U(l), so 



7~(-g.ym '■= L 2 (Asu(2), d[i S u(2)) ® L 2 (A u ^,d(j,u(i)) (38) 

In order to analyze the action of the Hamiltonian and to compute its kernel,it is 
convenient to introduce a complete orthonormal basis on the Hilbert-space (l38| ) . 
On the space of L 2 (Asu(2), d^su (2)) these are called spin network functions and defined 
as follows: let 7 G S be a graph and denote its edges and vertices respectively 
by (ei,...,ejv) and (vi,...,Vy). Associate a coloring to each edge defined by a 
set of irreducible representations (ji, • • • , jjv) of SU(2) (half-integers) and contractors 
(pi,...,Pv) to the vertices where pi is an intertwiner which maps from the tensor 
product of representations of the incoming edges at the vertex v\ to the tensor product 
of representations of the outgoing edges. A spin network state is then defined as 



N V 

\ T ( A )kip- <8>*(M^)) ■ <g)p* (39) 

i=l k=l 

where • stands for contracting at each vertex the upper indices of the matrices 
corresponding to all the incoming edges and the lower indices of the matrices assigned 
to the outgoing edges with all the indices of pu- 

In the case of L2(Au(i),dfiu(i)) one must simply replace SU(2) with U(l) in the above 
definition - these are called flux network functions ([9],[l0]). Since U(l) is a commutative 
group, we will have the following definition: for each edge of the graph associate an 
integer U. Then the flux network function is defined as 



\ F (A)) lX --=f[(K(A)) k (40) 
1=1 

What remains is to define the operators corresponding to the connection and the 
electric field on the Hilbert-space. If we want to implement the Poisson-brackets in the 
quantum theory in a diffeomorphism covariant way, we have to use smeared versions of 
these fields. In the case of gauge fields the natural candidates are the holonomy and the 
electric flux respectively: 
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h e (A) = Pexp A (41) 

J e 

E(S) = [ *E, (42) 



where e is a path and S is a surface in S, and *E is the dual of the electric field. 
Then the action of the corresponding operators will be defined as: 



k(A)f(A) := h e (A)f (43) 
E(S)f(A):=fh{E(S)J(A)} (44) 

3.2. Scalar field 

The crucial point of quantizing the scalar field is (see [5] or [H],[T2]) that the field 
should be real valued. In our case the original variables are complex, but this does not 
cause significant difficulties since we can introduce new fields which are real, thus the 
usual techniques can be applied on them. The only non-trivial problem is an additional 
ambiguity which arises because this can be done more than one way. What we are going 
to do is introduce two kinds of different choices for the configuration variables and the 
momentum operators. 

Case A: 

The most natural choice is to define the operators with the help of the real and 
imaginary parts of the fields. Let v be a vertex of a graph 7 with coordinates x v . Then 
let 



U(X,v) : = exp(zA3?($(x„))) (45) 
U(S,v) : = exp(i^($(^))), (46) 

where A and 5 are arbitrary real numbers which are required because otherwise the 
quantization would not be general enough (see pj] and [12J for details). The variables 
for the momentum operator should be (B is an open ball in E) 



11(B) = / d d x$t(n) 

J B 

U(B) = [ d 3 x%(7r) 

J B 

and thus the Poisson-brackets of the variables will be 
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{U(B),U(X, v)} = 5 vnB ^U(X, v) (47) 

iS - 

{U(B),U(5,v)} = -5 vnB , v -U(5,v) (48) 

{U(B),U(5, v)} = {ti(B), U(X, v)} = (49) 

The transformation rules of these quantities with respect to the (smeared) gauge 
transformation (G(A) = J d 3 xGA) are: 

{(7(A), U(X, v)} = ieAX%mx v ))U(X, v) 
{G(A),U(8,v)} = -ieA6$t($(x v ))U(6,v) 

Since we have two fields, the Hilbert space for the scalar field is a tensor product 
H sc = H(U)($$H(U), where the Hilbert spaces H(U) and H(U) are the linear 
combination of the following monomonials: let v = v±,...,Vn be the set of vertices 
for some 7 graph and let A and 6_ be two sets of real numbers, each pair associated to a 
vertex. Then a basic element of H(U) is constructed as follows: 



N 

\$ y = Y[U(\ k ,v k ) (50) 
fc=l 

In a similar fashion 



\^ = Ylu(Sk,v k ) (51) 

k=l 

will be a basic element is 7i(U). Both |A) 7 and |5) 7 form a complete orthonormal 
basis, that is y(A'|A) 7 = 8\x'S ltl ' an< ^ the same is true for |<5) 7 . 
Thus elements of 7i sc are linear combinations of of monomonials |A) 7 |5) 7 . 
The operators are defined in the same way as in the case of gauge fields: 



U(X,v)\X),:=U(X,v)\X), (52) 

n(£)|A) 7 := ih{U(B), |A) 7 } (53) 

U(6,v)\^ :=U(5,v)\l), (54) 

n(5)|I) 7 :=ih{T[(B),\l}^} (55) 
Because the U(l) group is commutative, the action of the operators are very simple: 

fl(S)|A) 7 = -|^A,|A) 7 (56) 

Vj eB 

#(A,v)|A) 7 = |A% 

K = K + S V ^X (57) 
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and similar expressions hold for U (5, v) and 11(B) . Also, because of (j4"9~|) 11(B) | A) 7 = 
fl(B)\l)^ = 0. 

Case B: 

Another way is to use the absolute value and the argument of $. Actually these 
are equal (up to constant factors) with the fields r] and 6 respectively, so we suggest 
the following operators for the multiplication operators: 

U v (\v) := exp(i\r)(x v )) (58) 

U @ (S,v) := exp(i5^^). (59) 
a 

For the momentum operators it is plausible to use the quantities 7r v and 7r instead 
of R(U) and Qf(II): 

Tl v (B) = / d 3 XTi v 
Jb 



Ile(-B) = a d 3 XTc e 
Jb 



The Poisson-brackets of these variables are a bit different then in case A: 

{U v (B),U v (\v)} =i^\5 vnB>v U v (\,v) (60) 

{U e (B),U e (5,v)} = i~85 vnB!V U e (8,v) (61) 
{U v (B),U e (S,v)} ={U e (B),U v (X,v)} = (62) 
The transformation rule for these variables with respect to gauge transformations 



arc: 



{G(A),U v (\,v)} = 

{G(A),U e (5,v)} = -UeA8&(x v )U e (5,v), 

which means that U v (X,v) is gauge invariant and the transformation rule for 
Uq(S,v) is 

U @ (5,v) » U @ (5,v)Ue(^,v)- 1 . (63) 

The construction of the phase space is completely identical to the construction in 
case A, the only difference is that one has to replace the old variables with the new ones. 
To avoid confusion, TL™ W = H(J7 ?7 ) <^)'H(Us) will stand for the new phase space, 



N 
k=l 



(64) 
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will label an element of 7i(U v ) and 



N 



\S @ ), = l[Ue(Sf,v k ) (65) 



k=l 



will be an element of "H(Uq). The action of these operators are completely the same 
as in case A: 



n fl (fl)|A") 7 = --5>;iA"> 7 (66) 

^(A,«)|A") 7 = |A rt? ) 7 

X? = X*> + S v , Vi X (67) 



5.5. Regularisation 

In order to quantize this system, one first has to rewrite the Hamiltonian in terms of the 
variables defined in the previous section - this is called the regularisation procedure. The 
key observation is ([5]) that if the gravitational field is dynamical, one can construct a 
well defined, diffeomorphism covariant Hamiltonian operator. In the article mentioned 
above the reader will find the detailed analysis of the gravitational, Yang-Mills, scalar 
and fermion fields. Since the method is quite lengthy, we are going to concentrate only 
on those terms that are different to the ones mentioned above. Specifically these are 
the terms that contain the scalar field. We will deal with the two kinds of description 
(the original case with $ and $* and the case with new variables r\ and 0) separately. 



Case A: 

Although later we will use the formulas involving rj and we shall provide 
the regularisation of the original Hamiltonian, since it has some non trivial steps. 
The potential term is the simplest: since $$* = 3?($) 2 + S($) 2 and 3?($) = 
arccos ( u (X> v ) + u — we can wr j^ e this (using the notations of [5J) in the following 
form: 



x 



, ■ x 

4 

v 

1 fl'iX.r) ■ U (A. f)\ 2 1 fU(5,v) + U- 1 (5,v)\' 



^ 2 arccos I I + — arccos I | - a 



(61 



One may wonder why we used the arccos function instead of e.g. the logarithm. 
The main reason is that since spontaneous symmetry breaking requires the ground state 
of the potential, we are forced to regularize the potential term to be self-adjoint. It is 
easy to see that the above operator is self-adjoint, but this would not be the case if we 
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used the logarithm function. Of course there are still ambiguities in the regularisation, 
but this certainly narrows down the possibilities. 

In a similar fashion, one replaces ttty* = ^R,(n) 2 + Q^tt) 2 in the kinetic term to obtain 

H P ^Ni V ) X ^ v ; W GM, (69) 

where X(v) and X(v) are the invariant vector fields on U(l) and Gi(v) contains 
only gravitational variables and it is the same as in [5]: 



r („\ - UK LMN v 

tri{v) — — — 9 T d f. > e e €ij k ei mn x 

v(A)=v(A')=v 

x Q l SLm (v, ~)Q? MiA <)(v, \)Q n SNm (v, ~), (70) 

where Q k e [v,r) = tr{rkh e [h~ l ,V(v) r ]), h e being the holonomy of the Ashtekar 
connection along edge e, V is the volume operator and Tfe are the generators of SU(2). 
The derivative term needs a more careful treatment. First we have to rewrite it in terms 
of and 



D a $(D b $)* = (d a + ieA a )(M(<5>) + - zeA fe )(9tf($) - z3($)) 

From this we can see that we need to regularize the expression (d a ± ieA a )^t(^). 
This is quite similar to the derivative term <9 a $ ± A a in [T], the only difference is that 
we have a iA a 3l(§) term instead of A a . Though this seems a minor change, it turns out 
that the regularized expression for this covariant derivative is more complicated, which 
is due to the fact that it contains the multiplication of the two fields. We can overcome 
this difficulty by doing the regularisation in a step-by-step way. First we note that for 
small At 

h s = l + ieAts a A a + o(At 2 ) 
for an edge s. This means that (v is the beginning of the edge s) 

(h s - 1) arccos ^ U ( X ^ V ) + ^(\v) ^ = ieX Ats a A a $t{$) + o(At 2 ), 
so if we take into account that 

U(X, s(At)) = 1 + iA(»($) + Ats a 8 a ^)) + o(At 2 ), 
we arrive to the following regularized expression: 
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U(X, s(At))[l + i(h s - 1) arccos ( 1 ' V) + ( ' V) )]E/(A, v)' 1 = 

= [1 + iA(ft($) + Ats a <9J?($))](l - eAAts a y4 a K($))[l - iA9fc($)] + o(At 2 ) = 
= 1 + *AAts a (<9J?($)) + 2eA a 3?($))) + o(At 2 ). (71) 

We obtain the same result for (d a + ieA a )^s(§) if we replace U with [7. Also the 
term (d a — ieA a )$t(&) is obtained by replacing h s with hj 1 . To simplify the result let 
us introduce a notation: 

W(v, s, A) = i[/(A, s(At))[l + - 1) arccos ^CM + ^CM ^ ][/(A, i;)- 1 
s, 5) = ^(<y, s(At))[l + i(h s - 1) arccos ^ ^M+^(M \ ^]u{5, v)' 1 - 
Using the fact that hj 1 = h s -i, the regulated expressions are the following: 

W(v, s, A) = iAts a {d a $l(<5>)) + ieA a $l(<S>))) + o(At 2 ) 
W(v, s-\ A) = tAts a (d a ^)) - ieA a M($))) + o(At 2 ) 
W(v, s, 5) = tAts a (d a %($)) + ieA a <Z(<f>))) + o(At 2 ) 
W(v, s-\ 5) = iAts a {d a %{<S>)) - ieA a %{<5>))) + o(At 2 ) 

This way the derivative term will be the limit of 

2 ^ E(v) 2 

x Yl [Wa(v, s n , A) + iW A (v, s n , S)][W A ,(v, s;\ A) - iW A ,(v, s" 1 , 5)]G™», 

t>(A)=u(A')=u 

where 

and the A and A' subscripts represent the tetrahedra where the holonomies and 
pointholonomies should be calculated. 

Case B: 

In this case one should be careful since the two scalar fields do not appear in a 
symmetric way. For instance, the potential term contains only rj, so one simply replaces 
77 = ± arccos ( MW^M j to obtgdn 
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Hp°t = \ N ( V )^^2 arCCOS ( 

v 

U v (X,v) + U-\X,v 



U v (\,v) + U- 1 (\,v) 



— arccos 
A 



1 2 



2a 



(72) 



The terms containing tt v and tvq can be treated in the same way as in the previous 



case, one just has to be careful since the later contains the expression 



(a+r)) 2 



But it is 



easy to see that if one carries out the regularisation procedure as in [5] the only difference 
will be a term which is the above fraction expressed with the variables U v . Thus the 
result for the two kinetic terms will be 



X v (v) 2 + 



i arccos f En^Mil^l , +a 



-2 



XM 2 



E(v) 



■om) 



The derivative terms for the two scalar fields are also different. In the case of r/, 
one only needs the expression (for small At): 

U V (X, s(At))U~\X, v)-l = i\Ats a d a r] + o(At 2 ), (74) 
thus it has contribution to the Hamiltonian 

v \ ' v(A)=v(A')=v 

x[[/,(A, Sr (A'))^- 1 (A,t;) - l]G™{v) (75) 

For the field G we remind the reader that in [1] the same kind of coupling appeared 
between the scalar field and the Maxwell field. Thus we may use the approximation 
mentioned there: 



U Q {5,s{At))h s U e {5 : v)- 1 - 1 = i5Ats\— 4 < A 



(76) 



Treating the term (a + ti) 2 as before, the regulated expression of this term will be 



#*r(e) = o£ 



N{v) 



1 f U v (X,v)+U- 1 (X,v)^ n2 

- arccos I ) + a 

-A V 2 / . 



v(A)=v(A')=v 



5 2 



[Ue&SniA^Ue^v)- 1 - 1] [U @ (5, s r (A'))h Sr U e (5, v)" 1 - l]Gf{v) 



(77) 
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In contrast to the Proca field, the mass here is represented by an operator, namely 

m(v) = e — arccos ( 1 + a . (78) 

-A V 2 J . 

Since we want to compare the two theories, it would be useful to work in a basis 
where U V (X, v) is diagonal and this is what we are going to do in this section. 

4-1. The spectrum ofU^X.v) 

Let |0) := |A^, . . . , AJJ ) be a base element for a 7 graph which has N vertices The 
action of U v (X,v) on this state is 

U V (X, v)\4) = | A* , . . . , X\ + A, . . . , K N )5(v, v k ) (79) 
This action suggests that we should look for eigenstates in the form 

|A"(A,,),A"):= ± (^M)lA'), (80) 

i— — 00 ^ ' ' 

where IX?) is an arbitrary state and A v (X,v) is a (yet) arbitrary number (this will 
be the eigenvalue for a given A at a vertex v). It is easy to verify that 



U v (X,v)\A«(X,v),>J>) =A(A >V )|A"(A >V ),A">. (81) 

We shall call these one vertex eigenstates because \A v (X,v), XP) is the eigenstate 
of only those U v (X,v') where v'=v.Since U v (X,v) is unitary, we can write A v (X,v) 
in the following form: A v (X,v) = exp(iA v (X,v)), where is real. In fact, since 
U v (0,v) = 1 and U v (Xi, v)[/ J) (A 2 , v) = U v (Xi + X 2 ,v), we obtain that A , '(A,f) is of 
the form A v (X,v) = T v (v)X. In summary, the spectrum of U v (X,v) is of the form 
exp(iXT v (v)), so instead of A v (X,v) we shall use T v (v). We can select an orthonormal 
basis from these eigenstates in the following way. Note that if there exists an integer n 
such that I A7> = U(X,v) n \Xl) then |P»(u),Aj> = e inXrv ^ {V 1 ! (v) , X>) . Because of this let 
us restrict ourselves to those (A' 7 ) that satisfy the condition < A^(A. Further more if 
we restrict the values of T' ? (f) so that < T v (v){^-, these states will form a complete 
orthonormal basis in the sense 



(rj( v ),A 1 |r2( v ),A 3 > = 

00 00 

= EE W> ■ ■ ■ > A ^ + kX > • • • ' A Ki A ^ ■ ■ ■ ' A "* 2 + a ■ ■ ■ » A S> >< 
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oo 
k=— oo 

= *(A^AJ)*(I>)-I>))), (82) 

where 5(AJ, Ajj) = 5(A^ 1)r/ - A^ 1)r? ) . . . ^A^ - A? "). To see that this is a complete 
orthonormal basis, one only has to check whether each original basis element can be 
expressed as the linear combination of the eigenstates. Let us suppose then that there 
exist complex numbers CxiT^^v)) such that 

I dT v (v)C\v(T v (v))\T v (v), \ v ) = |A"-> (83) 

A" J 

for each \X V ')- Because of orthogonality we obtain for the coefficients the following: 

oo 

C A ,(r») = (T T >{v),\ r '\\ r >') = (X v \ exp(-iifcAH(t;))^(A,7;) fc |A^ (84) 

k=— oo 

Now if for a \XP) there exists an integer n such that | A 77 ) = U(X,v) n \X v ') then the 
corresponding coefficient will be 

C7v>(r"(u)) = exp(-mAr»), 

otherwise it is zero. It is easy to see that this correspondence is unique and since the 
original basis is complete, we verified our statement. 

We define the graph eigenstate in a similar fashion. Let 7 be a graph and (A'') be 
an arbitrary state on that graph. For each vertex let T v (vi) i = 1 ... AT be a real number 
satisfying < P*^)^. Then the graph eigenstate will be the following: 



|r",A"):= ••• (e lfelAr ' ?( " l) ^(A,t;i) A;i )...(e i/cjvAr ' ,( ^ ) ^(A,t; 7 v) A;jv )|A' ? ) (85) 

k\=— 00 fejv=— 00 

Using the results obtained for the one vertex eigenstates we can find an orthonormal 
basis in the case of the graph eigenstates: if < AJJ(A for all v we get 

N 

(rimi%) = A3) II - rZM) ■ ■ • WK) - « (86) 

fc=i 

Also we can express any state in terms of graph eigenstates with the help of the 
following expression: 

[ dr>Cxv(r>)\r>, yi) = \x v >) (87) 

A" •* 
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where 



J dT 1 = J dT^{ Vl ) ...J dT>(v N ). 



What remains is the action of the momentum operators on an eigenstate. This is easy 
because of the following: 



k=— oo 



= (6„„X k + i6^,-^--)\r>(v'), A) (88) 
With a completely similar analysis one can show that 

X(v)\T\\) = (5 V)Vk \ k + A) (89) 

5. Solution to the constraints 

In [1] we sketched how one could solve the constraints of the theory. In this section we 
will follow the same procedures mentioned there - especially in the case of the scalar 
constraint. This method can also be used in this case, with one difference, namely 
that for the fields rj and we do not work in the usual basis, rather in the Fock- 
space. Since for the other fields the algorithm remains the same, we will concentrate 
only on the scalar fields. Solving the diffeomorphism- and gauge constraints will be 
rather simple, so we start with them. Then - in order to simplify things - we will 
introduce a compact notation where we separate the scalar fields from the others, which 
is described in Appendix A. This is motivated by the fact that the scalar constraint is 
quite complicated, but with the new notation the structure of the equation will be much 
easier to examine. 

Let us start with the diffeomorphism constraint. As it was pointed out in [5] the 
infinitesimal generator of the diffeomorphism constraint cannot be implemented in the 
quantum theory, thus the techniques used to solve the Gauss- or scalar constraint cannot 
be applied here. The strategy is to use group averaging to solve the constraint, which 
can be generalized to the case where matter fields also appear (see [5] for details). Since 
these are applied only to graphs not to the labels means that it is independent whether 
we use the Fock-space or the dust network space. 

The gravitational Gauss-constraint is the same as in [6], so we can solve it by restricting 
ourselves to gauge invariant spin network states. 

The U(l) Gauss-constraint contains variables of the electromagnetic field and the scalar 
field 6 so we analyze it in detail. The (smeared) integrated constraint 
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GA = / AiDaE. - aen e ) (90) 

J a 

can be regulated in the following way: Let us look for solutions in the form 

*= £ / Cc s ,/,a,5(L f) (s| (/| (r, A| (r, I| . (91) 

It can be verified that the quantum version of the above constraint is the following: 

(*l£ A »[£ / e -^ + i ^y)]l $ ) = ° ( 92 ) 

for all spin color network state |$). Here l e is the integer on the edge e (this comes 
from the flux network). Since A„ is arbitrary the above equation is equivalent to 

/ dr(vyC sJA! ?(L?)(r:i\ £ J.- (5 V + ^Jrr)\Ll) = 0, (93) 

J eDv=v ^ ' 

where we inserted (1911) to the constraint equation and used orthogonality of spin 
color network states. Now after partial integration we obtain a (functional) differential 
equation on the coefficients C s t x j(E, T): 

[£ le~(S v - ^—)]C S J M (T,^) = 0. (94) 

edv=v ^ ' 

Since we have a similar equation for all v, the solution to this constraint is: 

C s>fM (L,L) = C sjM (T) IJexp[-<( £ h- 5 v )T(v)], (95) 

v eHv=v 

where the coefficients C s f X g(T) are arbitrary. 
What remains is the scalar constraint. If we look at the Hamiltonian, it is clear that the 
constraint equation will be a differential equation with respect to the variable T. First 
we write down this equation. The condition we have to solve is 



{9\H\(f)) = (96) 

for arbitrary \(p). Again we can say that the support of N is at only one vertex v. 
Substituting (I9~TI) into the above equation we obtain 



s',/',A',I' 

x( s / |(/ , |(r',A'|(r / ,I / |F|r,I>ILA>|/)| s ) = o 
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In pQ we have shown a method (generalizing the results of [3]) which simplified the 
above equation by turning it into a finite number of equations. The main idea is that 
we take a basis element \s)\f) (discarding the scalar field for the moment) and we create 
a set containing basis elements appearing in H\s)\f). We continue this procedure 
and construct S^ n ' recursively from S^™ -1 -*. After this we search for solutions of the form 

(*i = E E c sM(fi 

i </|<s|eS» 

where s is a basis element. One can show that if we substitute this into the constraint 
equation we arrive to a finite number of conditions (details can be found in the mentioned 
articles). Now we apply these results to the gravitational and electromagnetic fields and 
we arrive to a finite system of linear differential equations. Since we concentrate only 
on the scalar field |T, A) at the moment and the Hamiltonian contains several terms, we 
shall calculate each term separately and introduce a compact notation. This notation 
is introduced in appendix A where the reader will also find the terms of the scalar 
constraint. The conclusion is that the scalar constraint is actually a system of linear 
differential equations of second order: 

E H ri{K - = - E HMDCAL) (97) 



v 



where 



Hn{T) = 

= L p(l*+$ + H " YM ^ + H ^(m + af + Hfi(v)r(vf(r(v) + 2a) 2 + 
+Hfj(v) exp(-2zAi») - Hf,j(v) exp(-zAr») + H?,j{v) (98; 



To simplify this term we look for solutions of the form 

Ci{v) ■.= c / (r) exp (-i n^W)' 



since 



The other terms will also contain a factor exp(— % Y[ v \ v T(v)) so this drops out of the 
differential equation, leaving us with the following formula for Cp(T) : 

E H "sTwA<r) = E HMDCAD (99) 



Spontaneous symmetry breaking in Loop Quantum Gravity 



21 



5.1. Solving the scalar constraint 

This system of linear differential equations can be solved using the method we shown in 
Appendix B if the matrix Hf,j is invertible. If it is not invertible then let us diagonalize 
the left hand side, i.e. find a unitary U such that UH P U _1 = diag(ki, . . . , fcjv) where 
fci,...,fcjv are the eigenvalues of H p . Let us order the eigenvalues in a way that 
hi, ... , ku (M(N) be all the zero eigenvalues. This means that the first M equation in 
this case is not a differential equation but only an algebraic equation. Since in this case 
the left hand side is zero, the right hand side is zero if and only Y2i>=i H'ri( v )Ci' = 
etc. for all matrices appearing in Hjij(T) (Hj?j(v) = (UH Q (v)U _1 )/'/ etc.), which 
means that after solving the algebraic equations we again arrive to a system of linear 
differential equations but with an invertible matrix on the left hand side. So from now 
on we consider Hf, T to be invertible. 

To have a correct solution we must specify the initial condition on C/'(r) and g^Cr(r). 
The fact that T(v) = —a can be interpreted as the disappearance of the field rj implies 
that 



E tf/'/(r)Cj'(r)]|r ( t,)=-« = o (ioo) 



i 



be the first condition. With the same reasoning the second condition is that the 
momentum of the field should disappear. In this case (since all A„ are zero) we arrive 
to the condition 



[^CAT)] lr(v)= _ a = (101) 

Since Hn has a complicated structure, the differential equation cannot be solved 
explicitly. However we can solve it in some special case. 

First let us consider the case when T(v) —a. In this case the system of differential 
equations takes the form 

Now if we multiply both sides with (H p ) _1 and define bj = ((H p )- 1 H G+YM C(-a)) / 



we get 



Ci{T)- — \ V , 9 C f (T) = 6 J (103) 



5T(v) 2 JW (I»+a) 2 
The general solution of this differential equation is the following: 



dliL) = { l- ] LM2 bl + CKT{V) + ar + CKnv) + a) " 2 ' (1 ° 4) 
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where 



x _ 1± y/l + 4L(v) 2 




and C),Cf are arbitrary constants. From L(v) 2 > follows that n\ > 1 and rz 2 < 0, 
which means that if L(v) ^ then C/QT) is singular in T(v) = —a. 
Now let us consider the initial conditions. If L(v) = then 



Cj(T) = (r( ^ 2 +Q)2 & z + C£(r(v) + a) + C 2 . (105) 

Substituting into (11001) and (110 II) implies that C] — and 6/ = 0. Further more 
from the definition of bj comes that bj = ((H P ) _1 H +rJW C 2 )j, so the solution is: 



CAT) = Cl (106) 
where C 2 must satisfy the condition 



Hpf YM Cj = 0. (107) 

This is not a surprising result since if L(v) = then substituting this into the 
constraints we obtain a theory completely equivalent to the electromagnetic field coupled 
to gravity. If we rewrite the scalar constraint of this theory in terms of the notation 
used in appendix A, we obtain the above condition. 

What happens if L(v) ^ 0. In this case C 2 = so that the solution does not become 
singular at V = —a. Substituting into (110 II) will yield the identity 0=0, so we must 
check (II 001) . For L(v) = ±1 this will be singular so in this case C\ = and only the 
first term survives, but it will be zero too. Thus in this case the solution near F = —a 
is zero in first order. For |L(u)|)l the condition fllOOj) is also an identity. But in this 
case Cj (—a) = 0, so bj = 0, thus 



C/(E) = Cj{T{v)+a) ni . (108) 

This solution tends rapidly to zero as T — >• —a (especially if L(v) is large), so as 
we reach this limit, the amplitude of the solution coming from the L(v) = case will 
become significantly larger. In fact the larger L(v) is, the amplitude becomes much 
smaller in this region. So we can say that if T(v) + a « (which - as we will see later - 
can be interpreted as the mass is about zero) states which for which L(v) = have the 
highest probability while the larger the smaller this probability will get. 

These results show that in contrast to the Proca field, this theory provides us with the 
different amplitudes for different masses. However because the two theories are - in some 
aspect - very similar, it would be desirable to provide the solutions of this theory which 
can be identified as the solutions to the Proca field. The basic idea is very simple: we 
compare the two Hamiltonians. If we look at the matrix (1981) in our differential equation, 
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in the case T = it will be the same as the Hamiltonian of the Proca field. So one just 
needs to imply the conditions 

A'(E) = (109) 



5T(v) 2 

E#/'/(r)cv(r))ir=o = o (no) 

I' 

The problem is that in this theory this will provide a distributional solution in the 
following sense. In the case of the Proca field the mass is fixed, which means that we are 
interested in solutions where T is constant. But now we have a differential equation so 
T is continuous. The way out of this is we say that in the interval (— e, e) we solve (I109p . 
and outside this interval CV(r) is zero. The required solution will the limit e — * 0. The 
reason for this strange behavior is that the equation we gained looks not like the Proca, 
but the linear combination of all the Procas. 

6. Mass 

In quantum field theory the mass is the coefficient of the term in the Hamiltonian which 
is quadratic in the boson field. However in this case we shall define the mass as an 
operator corresponding to the classical expression n + a. The reasons for us to do so 
are the following: First - as was shown at the end of section 2.2 - the expression r] + a 
corresponds exactly to the mass parameter of the Proca field (The term (77 + a) 2 not 
only appears in front of the quadratic term of the bosonic field but also appears in the 
denometer of the kinetic term of the other scalar field). The second reason is that in 
this case we can simplify our analysis regarding the scalar-bozon interaction. This new 
interpretation - as we will see - gives a better understanding of the mass generation 
in the Hamiltonian framework. Note also that the substitution 77 = gives back the 
"original" mass. 

Let \ty) := Y^x v I diyCxv^Yy)]!^ 1 ,)^ 1 ) be a solution of the constraints. Then we can 
define the "mass operator" as 

m W) = — arccos ■ M/> = 

A V 2 /' 1 

= Y1 [ C^(r 7 )r>)|r 7 ,A ,? ). (111) 

It is clear from the definition that this operator is self adjoint, thus it has real 
eigenvalues. Further more its spectrum is continuous. Its expectation value is 

m {%v) = W c^q,(r?)c A ,(rl) e (r^( u ) + a)(r? ! A?ir^,Ai) = 
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= W dr\C^(D\ 2 e(^(v) + a), (112) 

thus for a graph 7 we may define the mass as 

m(¥,7) = X>CM = X) / Cl^(H)| 2 ^e 2 (r>) + a). (113) 

i> A* 1 u 

This means that r(v) + a can be interpreted as "mass in a vertex". If we look at 
a state where all F are zero - the "vacuum" (note that it is NOT the usual vacuum 
since we are not in the Fock space representation) - we obtain states with mass ea. 
But one may ask whether this is an observable or not. If one checks the commutator 
of the constraints and rh the only non- vanishing term will be the [Hp, rh] commutator, 
which is proportional to X(v). This means that if take the subset of the solutions where 
X(v)ty = 0, the mass operator will be an observable. But if we look at the action 
of X(v) in our new basis in (IBB]) one will find that this is equivalent to the condition 
(11091) . So rh is an observable if the solutions are those which are equivalent to the 
solutions of the Proca field. But one may say that there are other solutions as well, 
since one does not have to impose (11101) . The answer is that these states are special 
cases which are contained in the Proca solution. This is because in this case one has to 
solve Yip Hi'i(T)Ci'(T) for all T, which means that the solution will have to be in the 
kernel of all matrices appearing in Hpj. 

All in all the mass operator is an observable if if the solutions are those which are 
equivalent to the solutions of the Proca field. Since the Proca field did not have a 
potential term, the correspondence is correct only if we consider states where all T are 
zero (note that other states the mass operator is also an observable, but it describes 
interactions). 

7. Summary and open questions 

In this paper we analyzed the mass generation to a U(l) vector field via spontaneous 
symmetry breaking in LQG and compared the results obtained for the Proca field. Even 
at the classical level - after introducing new variables 77 and - the two theories had 
many similarities. The main difference was the extra scalar field and the potential term 
in the case of the spontaneous symmetry breaking, and where the Proca field had a mass 
parameter, we obtained a field. Thus it was not a surprise that the quantized theories 
were also similar, and in the case of spontaneous symmetry breaking mass became an 
operator. We defined a new basis in the quantum region, where the motivation was 
to find the eigenstates of the configuration variable of the scalar field. By choosing 
this new basis we were also able to rewrite the constraints to finite linear systems of 
differential equations, thus we were able to analyze only the scalar field dependence of 
the theory. We were able to (partially) solve the constraints and describe the behavior 
of the states when T + a (i.e. the mass) tends to zero. We found that there exist states 



Spontaneous symmetry breaking in Loop Quantum Gravity 



25 



which are non-degenerate in this region, and further more at T + a = there is only one 
non-zero amplitude, the one which belongs to a solution to the "gravity coupled to the 
electromagnetic field" case. 

The eigenvalues of the mass operator are continuous (though they have a discrete 
structure due to the discreteness appearing in the coefficient matrices) and real, but not 
necessarily positive (one needs extra input for this). A very interesting result is that 
the mass operator is an observable if the states are in the kernel of the corresponding 
momentum operator. This extra condition implies that in this case the other scalar field 
also becomes a gauge. This means that if we want a physically relevant mass operator, 
the scalar fields will not be real particles. 

In the light of our results we can claim that - though the two theories are very similar - 
the spontaneous symmetry breaking has more advantages: 1. We are able to calculate 
the mass dependence of the states without solving the entire theory 2. We can produce 
the limit m — > without difficulty (the states have non-singular solutions) 3. In the 
case of the Proca field for different m we have different theories, while the spontaneous 
symmetry braking deals with all values. This is important if we want to calculate 
transition amplitudes between states that have different masses. 4. The mass is an 
eigenvalue of an operator which can be an observable, while in the case of the Proca 
field mass is a parameter. 

In our analysis it was crucial that the scalar field had a commutative group, otherwise 
the eigenvalues of the configuration operator would have been hard to find. Thus it is an 
interesting question that in the case of non-commutative groups how can we generalize 
these results? But the commutative case also provides a few questions, like the complete 
analysis of the differential equation (l98l) . The main question is what kind of restrictions 
do we have to make on the coefficient matrices to have a well defined, non-singular 
square-integrable solution. 
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9. Appendix A 

Here we introduce a notation which simplifies the scalar constraint. Let us 
introduce a multi-index I for the indices s,f, F,5 so that a type of expression 
(s'\(f'\{T,6\X\T,6_)\f)\s) will be denoted as Xn. Now consider those terms that do not 
contain U v or X v . These are H gra v and Hym, the Hamilton operator of the gravitational 
and Maxwell field. So in our new notation the contribution of these terms to the 
constraint equations will be the following: 
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s'J'A'A' 



/ dr(v)' / dr(v)XH^ ff ,5^5 M 5(r(v)', T(v))S(T(v)', T(v)) + 



+^ / 7 ; GL^%^(^(^;)^^(^))5(^(^;)^^(^;)))^, / , A ,,K^^^') = 

= X)flff y »^(r), (114) 

/' 

where we performed the integration and sum on the variables related to the scalar 
fields. 

The terms containing U v or X v will be treated as follows: we shall write the dependence 
of these fields explicitly, while other expressions will be denoted (using the short 
notation) as 0}, T etc. For example the contribution of the potential term H pot will 
be denoted as follows: 



J2 J dv{v)' J dr^'ijv^ii^l^laO^^f.n^^l 

s'J'A'A' 

xC s>j>,\',5'(T >H)<r, A| arccos ^ ^— £ J x 



x 



x[i arccos (^M±^(M) + 2 a]V,A'> = 



■ u v (\,v) + u- l (\,vy 

2 

= £ ^>)CHr)I>) 2 (I>) + 2a) 2 (115) 

7' 

The derivative term Hd e r(®) contains the operator U v , so we have 



./ // \ f & j 



x 



2£M 2 

1 /C/ 7? (A^) + C/- 1 (A,^x 



x 



(T, A| [i arccos (W> + »« ^) + „] |r', A') X 



1 .= 



x X) ^(f,l|[^e(5,^(A))/i s ^e(5^)- 1 -l]x 

d(A)=d(A')=d 

x[[/ ( ( 5, Sr (A'))/, Sr -7e(5^)- 1 - l]|r', ^Xsl^Cr;)^')^, y,r(I', l') = 

= ^^T 1 Cp(r)(r(^) + a) 2 (lie) 

/' 

The other derivative term is a bit trickier since the U v is evaluated in different 
vertices. We have 



i dV ^l df ^y\§^ E ^(r,A|[K ; (A, Sn (A))^ 1 (A,t;) 

s',f',X',S' { ' v(A)=v(A')=v 

x[u v (\,s r (A'))u-\x,v) - i]|r,A')( s |GTMk')^, r ,y,r(r , 'r / ) = 
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= J2( H ri e M- 2iXr ( v )) ~ ffj?jexp(-iAI») + H^^C^r), (117) 
v 

where the values exp(zAr(s r (A'))) etc. have been assimilated in the coefficients 
Hpj etc. since our differential equation will depend only variables in vertex v. 
The last contribution is the momentum term. It is convenient to use equation fl94|) so 
that one can simplify this expression in the following way: 



jdT{v)'jdt{ 

i f \> x' 



1 f U v (A,v) + U- 1 (\,v)s 1-2 



vy^^mx\x r ,(vr\r\^ J! 5(r(v)\r(v)) + 



s',f',X,S 

+ (r,A|^arccos(^ ~" v 1 ~ ' — ' r > — ~ ' ) + a\ )|r',A'} x 
x(fi|Xe(^) 2 )|f^l / )( S |G 1 (^;)| S ')5 // ,C' s , / , A , i K^',^ / ) = 

s',f>,\',s! {V> 

where L(v) = J2 V b is the sum 
10. Appendix B 

Here we describe the method to solve a system of differential equation of the form 

5(t) = H{t)c{t), (119) 

where H(t) is a N x iV matrix. 
The method is similar to the one used in cases of path ordered integration. First we 
integrate the equation: 

'c{t)= [ dhU^cXh) + c 1; (120) 



where 

cx = a;o). 

Another integration yields 

c(t)= [ dh [ dt 2 n{t 2 )c{t 2 ) +cit + c , (121) 
Jo Jo 

where 

c = c(0). 

Now we iterate this equation and arrive to the result 
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oo rt 



dt 2 ... 



+{t+Yl / dti / dt2 



00 i-t 
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